Direct measurement of spatial Wigner function with area-integrated detection 
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We demonstrate experimentally a novel technique for characterizing transverse spatial coherence using the 
Wigner distribution function. The presented method is based on measuring interference between a pair of rotated 
and displaced replicas of the input beam with an area-integrating detector, and it can be superior in regimes 
when array detectors are not available. We analyze the quantum optical picture of the presented measurement 
for single-photon signals and discuss possible applications in quantum information processing. 
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The problem of characterizing spatial coherence of opti- 
cal fields appears in a variety of applications. In biomedical 
optics, spatial coherence is being studied as a means to re- 
veal the internal structure of scattering media, thus providing 
novel modes of tomographic imaging. 1 1] In optical commu- 
nications, spatial coherence of light sources is necessary to 
achieve efficient coupling into fiber systems. This aspect be- 
comes even more critical in quantum communications, [i^SBl 
where non-classical sources of radiation are used, and detec- 
tion is often based on subtle interference effects that require 
good mode matching. 

In this paper we present a novel technique for characteriz- 
ingspatial coherence of light in the Wigner representation. 1 5 , 
00] Our technique has several features that make it distinct 
from previously proposed schemes. First, it is based on a se- 
ries of measurements performed with an area integrating de- 
tector that collects all available light power, rather than on 
spatially resolved detection employed is some of the previous 
proposals. 1 8, 9] Therefore, it can be used in regimes where 
array detectors are not available, for example in mid- and far 
infrared, or when single photon signals are involved. In such 
cases, our technique provides the optimal signal-to-noise ra- 
tio, analogously to the case of Fourier transform spectroscopy 
which is more efficient than a spectrometer with a single scan- 
ning detector. Compared to the recently proposed hetero- 
dyne schemes for measuring the Wigner function. iflCt Hill our 
method is self-referencing and it does not require any auxil- 
iary sources of radiation. 

Our technique is based on the following observation. For 
a complex quasi-monochromatic field E(£) varying in the 
transverse plane parameterized with the two-dimensional vec- 
tor the Wigner function is defined according to 00] 

W{x, k) = -L j d 2 Ze 2ik *{E*(x - £)E{x + £)>• (1) 

Here k is the transverse spatial frequency vector and the angu- 
lar brackets denote a statistical average. The above definition 
can be easily transformed to the form: 

W(x, k) = ^(^J d 2 i[e-^E{x - £)]*e ihi E{x + $)J . 

(2) 



This expression shows that the Wigner function at a point 
(cc, fc) is given by a spatially integrated overlap of the field 
e lk ^E(x + £) with a complex conjugate of its replica rotated 
by 180°, which corresponds to the transformation £ — > — £. 
The field e *E(x + £) can be obtained from the original in- 
put by two simple transformation: displacement in space by 
x and changing the direction of propagation by k, which in 
the paraxial approximation results in multiplying the field by 
a phase factor e lk ^ . 

The above idea can be implemented to measure the Wigner 
function in a setup depicted in Fig. [Da). The input beam is 
steered using the mirror Ml into a three-mirror Sagnac inter- 
ferometer. The displacement and the tilt of the mirror control 
the point at which the Wigner function is measured. The 50:50 
beam splitter BS generates a pair of replicas of the displaced 
and tilted input beam travelling in the opposite directions. The 
interferometer contains a Dove prism whose base forms 45° 
with the plane of the interferometer. Transmission through 
the Dove prism rotates each of the counterpropagating beams 
by 90°, which adds up to the 180° relative rotation required 
to measure the Wigner function. The operation of the setup 
can be most easily understood by analyzing propagation of 
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FIG. 1 : Experimental apparatus for measuring the Wigner function. 
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the replicas of an exemplary image through the interferome- 
ter, as shown in Fig.^b). The counterpropagating replicas of 
the input field are recombined at the output port of the inter- 
ferometer, and the emerging beam is focused with the help of 
a lens on an area-integrating detector. 

The intensity I recorded by the detector can be decomposed 
into a sum of three terms, I = h + h + Ii2- The first two of 
them are given by the space-integrated input intensity \E(£) 2 
and they remain constant as long as the aperture of the inter- 
ferometer does not clip off any of the input field. The third 
term originating from the interference between the coun- 
terpropagating replicas is proportional to the right hand side 
ofEq. 0: 

h 2 = ^^W( X ,k) (3) 

thus giving the value of the Wigner function. Here <p is the 
relative phase between the interfering replicas. The complete 
Wigner function can be therefore scanned by measuring the 
detector photocurrent as a function of the position and the tilt 
of the steering mirror, and subtracting the constant pedestal. 
We note that the factor relating the interference term to the 
Wigner function depends on the relative phase cp between the 
interfering fields, and it can assume in principle either a posi- 
tive or negative sign. We also note that if we extract the beam 
leaving through the second port of the interferometer, for ex- 
ample with the help of an optical circulator, then the constant 
pedestal can be removed by subtracting photocurrents mea- 
sured at both the output ports of the interferometer. 

We used the presented technique to measure the Wigner 
function for several beams derived from a He-Ne laser. For 
simplicity, we restrict ourselves only to one spatial dimen- 
sion, parallel to the plane of the interferometer. The trans- 
verse distributions of the electric field for the prepared beams 
were factorizable in the directions parallel and perpendicular 
to the interferometer plane, so that the input field can be effec- 
tively considered as one dimensional. The steering mirror was 
mounted on a pair of translation and rotary stages (Physik In- 
strumente, models M-014.D0 and M-036.D0) controlled by 
a computer, which also recorded the value of the photodiode 
current. The interferometer was constructed using 1 inch di- 
ameter mirrors and a beam splitter, and a Dove prism 15 mm 
in cross-section. The beam emerging from the interferometer 
was focused with an / = 25 mm lens on a detector with active 
area 3.6 mm x 3.6 mm. 

In Fig. |2] we depict the measured Wigner functions. The 
first two plots were obtained for a laser beam filtered through 
a single-mode fiber and (a) collimated or (b) made lightly di- 
vergent. These data show Gaussian Wigner functions with 
diffraction-limited widths along the principal axes. The tilt 
of principle axis in Fig.|2jb) reflects positive correlations be- 
tween the position and the direction of propagation character- 
istic for a diverging beam. Fig.|2jc) shows the Wigner func- 
tion for an input obtained by a reflection from the front and the 
back surfaces of glass wedge. The measured Wigner distribu- 
tion function shows two slightly converging Gaussian beams 
with different central transverse wave vectors. The oscillat- 
ing pattern located between the Gaussians is a signature of 



their mutual coherence. We fitted our experimental data to a 
theoretical model of a Wigner function describing a pair of 
partially coherent Gaussian beams, which yielded the degree 
of coherence equal to 1.02 ± 0.04. An additional feature seen 
in Fig. |2jc) is the presence of ripples on one of the Gaussian 
peaks, at twice the spatial frequency of the fringe pattern at 
the center of the distribution. These ripples can be attributed 
to interference with a third, much weaker beam generated by 
a secondary reflection from the glass wedge. 

In the theoretical analysis presented above, we assumed that 
the range of transverse spatial frequencies in the input beam 
is narrow enough to make diffraction effects across the setup 
negligible. In a geneneral case, the measured Wigner func- 
tion describes the distribution of the electric field in the plane 
in which the two replicas experience relative rotation. It is 
therefore critical that the replicas counterpropagating in the 
Sagnac interferomer are inverted at the same distance from 
the beam splitter BS that generated them. Diffraction after the 
Dove prism is irrelevant, as it does not change the total inten- 
sity which is the only quantity measured on the output. Also, 
when diffraction cannot be neglected, the parameterization of 
the Wigner distribution function obtained from the displace- 
ment and the tilt of the steering mirror needs to be compen- 
sated for the propagation distance from the steering mirror to 
the inversion plane. The compensation is given by a simple 
linear transformation of the Wigner function variables. |6j|7|] 

One important regime of optical measurements where ar- 
ray detectors are not easily available are single-photon sig- 
nals. When the output of the interferometer is measured us- 
ing a detector with single-photon sensitivity, for example an 
avalanche photodiode, the presented setup can be used to 
characterize transverse coherence of the output produced by 
single photon sources. J2J |3J] Development of such sources 
is presently a subject of extensive research effort owing to 
prospective applications in quantum cryptography. As se- 
cure quantum communication is most likely to be based on 
single-mode fiber links, efficient matching of the modal struc- 
ture of the source to the fiber is necessary to achieve high 
transmission rates and close options for eavesdropping at- 
tacks. 

The Wigner distribution function defined in Eq. Q charac- 
terizes second-order statistical properties of fluctuating clas- 
sical electromagnetic fields. However, it admits also a quan- 
tum mechanical interpretation which has a number of inter- 
esting consequences for our work. If the input signal con- 
sists of a single quasimonochromatic photon, then E(£), when 
properly normalized, is the wave function of that photon 
for the transverse spatial degrees of freedom. fl2l This cor- 
respondence has been recently used to propose encoding of 
information into spatial modes of a single photon character- 
ized by different orbital angular momenta. 1131 Separation of 
these modes was performed by interference between mutu- 
ally rotated replicas in a manner similar to that underlying our 
scheme. 

Following the quantum mechanical picture, the spatial de- 
gree of freedom of a photon is another optical realisation 
of a continuous-variable quantum system, alternative to the 
quadrature of a single light mode. There exists an analogy be- 



3 




FIG. 2: Surface and countour plots of Wigner distribution functions measured for (a) collimated and (b) slightly divergent single mode beams, 
and (c) a laser beam reflected off a glass wedge. 



tween the scheme proposed here and a direct measurement of 
the quadrature Wigner function of a light mode. ll4lll5ll In the 
latter scheme, the photon number parity operator is the equiv- 
alent of the spatial rotation between the interfering replicas 
that yields the Wigner function in the present case. This anal- 
ogy has interesting consequences when generalized to com- 
posite systems. Consider a pair of photons with perfectly cor- 
related positions, described by a pure entangled wave func- 
tion. Such a pair exemplifies the Einstein-Podolsky-Rosen 
paradox, and it forms a source of entanglement. For exam- 
ple, it can be used to test Bell's inequalities by sending the 
two photons to a pair of Sagnac interferometers and measur- 
ing the joint spatial Wigner function at selected points of the 



phase space. Such a measurement reveals quantum nonlocal- 
ity, analogously to the joint measurement of displaced parity 
operators on a pair of quadrature entangled beams that has 
been discussed previously.! 16] 

Concluding, we presented a method for measuring trans- 
verse spatial coherence of optical beams in the Wigner repre- 
sentation. The scheme is particularly beneficial in situations 
when array detectors are not available. It has a number of 
implications in quantum optics and quantum information pro- 
cessing: it can be used as a diagnostic tool for novel sources 
of nonclassical radiation, and also probe spatial entanglement 
in optical systems. This research was supported by ARDA 
Quantum Computing Program. 
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